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Abstract 

We study an anisotropic inflation model with a gauge kinetic function for a non- 
abelian gauge field. We find that, in contrast to abelian models, the anisotropy can 
be either a prolate or an oblate type, which could lead to a different prediction from 
abelian models for the statistical anisotropy in the power spectrum of cosmological 
fluctuations. During a reheating phase, we find chaotic behaviour of the non-abelian 
gauge field which is caused by the nonlinear self-coupling of the gauge field. We 
compute a Lyapunov exponent of the chaos which turns out to be uncorrelated with 
the anisotropy. 
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1 Introduction 



In an inflationary scenario, quantum vacuum fluctuations during inflation accounts for 
the origin of the large scale structure of the universe. The nature of such primordial 
fluctuations is understood from the approximate symmetry in quasi-deSitter inflation. First 
of all, in order to have the inflation, we need an approximate translation invariance for an 
inflaton field, which forbids strong nonlinearity in the action. Hence, we have Gaussian 
statistics of fluctuations. Moreover, the approximate deSitter symmetry makes the power 
spectrum of fluctuations scale invariant and statistically isotropic. Note that the scale 
invariance originates from the temporal part of the deSitter symmetry and the statistical 
isotropy comes from the spatial part of that. Since these predictions are based on the 
symmetry, they are robust in the standard single field inflationary scenario. 

However, the symmetry in quasi-deSitter inflation is not accurate from the point of 
view of precision cosmology. In fact, if the translation symmetry is exact, inflation never 
ends. Since the violation of the translational symmetry is characterized by the slow roll 
parameters, it is natural that the non-Gaussianity in a single inflaton model is of the order 
of the slow roll parameters ||T|. Apparently, the violation of deSitter symmetry is also 
characterized by slow roll parameters of the order of a few percent. Due to the violation 
of temporal part of deSitter symmetry, a deviation from the scale invariant spectrum can 
be expected to be of the order of the slow roll parameters. Actually, this deviation has 
been observationally confirmed B. In this line of thought, it is legitimate to suspect 
that the spatial part of deSitter symmetry also breaks down slightly. The violation of the 
spatial part of deSitter symmetry would lead to an anisotropy in the cosmic expansion, 
namely, anisotropic inflation. As a consequence, quantum fluctuations generated during 
the anisotropic inflation must have the statistical anisotropy. Therefore, it is quite natural 
to expect the statistical anisotropy of the order of the slow roll parameter. 

Historically, there have been many attempts to construct anisotropic inflationary mod- 
els p|-|T0f . However, it has been shown that these models suffer from the instability |jll|| , or a 
fine tuning problem, or a naturalness problem. Recently, a successful anisotropic inflation- 
ary model has been proposed More precisely, it turned out that the presence of a 
non-trivial gauge kinetic function in supergravity theory can accommodate an anisotropic 
inflation. It is well known that the supergravity is characterized by a superpotential, a 
Kahler potential, and a gauge kinetic function. These functions should be constrained by 
comparing predictions of inflation with cosmological observations. For example, the tilt 
of the power spectrum has given interesting information for the superpotential and the 
Kahler potential. The information provides a hint to the fundamental theory. Amazingly, 
so far, the gauge kinetic function in supergravity has been neglected in making predictions 
of inflation. The reason is partially due to the cosmic no-hair theorem which states that 



the anisotropy, curvature, and any matter will vanish once the inflation commences [ 14| . It 



has been proved that this is merely a prejudice [T^,|T3|- The reason is simply that we do 
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not have a cosmological constant because of the violation of the translation symmetry for 
the inflaton. In spite of the absence of the cosmic no-hair theorem, since the anisotropic 
inflation is an attractor, the predictability of the model still remains JIB-Hi]. Indeed, the 
imprints of the anisotropic expansion could be seen in the CMB Jl9| . 



From the particle physics point of view, it is important to explore the role of the gauge ki- 



netic function in inflation. In a previous work, we have considered an abelian gauge field [|12 
However, in reality, we have non-abelian gauge fields in particle physics models. Hence, 
the main purpose of this paper is to investigate a cosmological role of non-abelian gauge 
fields in an inflationary scenario. There are two important differences between abelian and 
non-abelian gauge fields, that is, the non-abelian gauge fields have multi-gauge-components 
and nonlinear self-couplings. Thus, in this paper, we focus on consequences stemming from 
these two features. 

Here, we should note other works on the statistical anisotropy. The statistical anisotropy 
generated by vector fields is first investigated using 6N formalism in [|(| where the possi- 
bility that the anisotropy appears strongly only in the non-gaussianity is pointed out. The 



model has been further extended in various ways [21-25]. In particular, the formalism has 



been generalized to non-abelian gauge models [26-28]. As a different approach, there are 
attempts to see the remnant of the universe before inflation ][29|-|34|| . This could be possible 
if the duration of inflation is sufficiently short. In the non-inflationary scenario, there is 
another mechanism for producing the statistical anisotropy |^SJ . 

The organization of the paper is as follows. In section II, we introduce inflationary 
models inspired by supergravity where the Yang-Mills field couples with an inflaton. In 
section III, we first solve basic equations numerically and obtain solutions which exhibit 
anisotropic expansion and chaos during reheating. Next, we present analytical formula 
for the degree of the anisotropy of the cosmic expansion during inflation. We also discuss 
observational implication of our finding. In section IV, we calculate a Lyapunov exponent of 
the chaos during reheating and find no correlation between the anisotropy and the Lyapunov 
exponent. The final section is devoted to conclusion. In the appendix A, we explain how 
to reduce the degree of freedom of the non-abelian gauge fields using the symmetry in the 
system. 



2 Inflation model in supergravity 

In this section, we present an inflationary model based on supergravity where we have 
a non-trivial gauge kinetic function for a gauge field. Although the gauge group could 
be general, we choose SU(2) for concreteness. Using Pauli matrices a a , we can define 
generators of SU(2) by T a = a a /2 (a = 1, 2, 3) satisfying 

[T a ,T b ] = ie abc T c , ti(T a T h ) = ^5 ab , (2.1) 
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where e abc is a Levi-Civita symbol and tr denotes the trace of the matrix representation. 
Here, 5 ab is a usual Kronecker delta. The 577(2) gauge field is defined as A = A^dx^ = 
A a ^T a dx^. We note that the gauge field A = A a T a has multi-gauge-component A a . 

The action for the gravitational field, the inflaton 0, and the gauge field reads 



S — d A xyf^g 



^- 2 R ~ \{dct>? ~ V{<j>) - l -f\<P)tr{F, v Fn 



(2.2) 



where R is the scalar curvature, g represents a determinant of the spacetime metric, V(4>) 
is a potential for the inflaton and the field strength F^ u of the S77(2)-gauge field is defined 
as F^ u = d^Ay — dyA^ + igylA^, A v ]. Here, gy is a Yang-Mills coupling constant. The above 
action is invariant under the local SU (2) gauge transformation, 

A, 7 -% 7 - —7-^7 , (2.3) 
9y 

where 7 G SU{2). The gauge kinetic function /(0) will be specified later. Equations of 
motion derived from the action Q2.2p are given by 



= V M 0V,0 - i^(V0) 2 - g^V(cj>) + 2/ 2 (0) tr[F w F/ - i^F 2 ] , (2.4) 

V 2 0-n0)-/(0)/'(0)tr(F 2 ) = O, (2.5) 
7J,[/ 2 (0)FH = , (2.6) 

where G^ u is the Einstein tensor, V M represents a covariant derivative with respect to the 
metric g^ u and we have defined the derivative ' = d/d<fi and the gauge covariant derivative 
D„ = V M + i^K, * ]. 

Now, let us consider a cosmological background spacetime. For simplicity, we consider 
the axially symmetric Bianchi type-I metric 

ds 2 = -dt 2 + e 2a ^[e-^dx 2 + e 2 ^\dy 2 + dz 2 )} , (2.7) 

where a describes the average expansion and a characterizes the anisotropy of the expan- 
sion. The symmetry in this spacetime is characterized by Killing vectors d x , d y , d z and 
£0 = — zd y + yd z . In particular, ^ generates the rotational symmetry in (y, z)-plane. Im- 
posing the symmetry on the inflaton and the gauge field A, we can reduce variables into 
the following form 

= 0(t) , A{x?) = vxit^dx + v 2 (t)(T 2 dy + T 3 dz) . (2.8) 

The gauge field A is parametrized by two functions, v i(t) and v 2 (t). In the appendix [A], we 
explain how to achieve the above form for the gauge field by using the symmetry. Substi- 
tuting Eqs. (|2.7|) and ( |2.8|) into Eqs. (|2.4| - |2.6|) , we obtain basic equations for the cosmological 
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background spacetime. From the time-time component of Einstein equations, we obtain a 
constraint equation 



3 , . 2 .i\ 1 ; 
1 



(_a 2 + < 7 2 ) + -0 2 + y(0) 



2 / 2 (e" 2a+4CT ^ + 2e- 2a - 2a v 2 2 + 2g^e- ia+2a vy 2 + c^e" 40 " 4 ^ 4 ) = , (2.9) 

where we defined a derivative ■ = d/dt with respect to the cosmic time. We also have the 
evolution equations 

-^a + %(a 2 + a 2 ) + ]-4> 2 - V 

K z K z 2 

+ ^f 2 (e- 2M v 1 2 + 2e~ 2a - 2 °v 2 2 + 2g 2 Y e- 4a+2 ° 'v\v\ + ^e" 40 " 4 ^ 4 ) = 0(2.10) 



and 



H 2 



Z. i QA.J- '" -2ct+4cr • 2 — 2a-2a • 2 „2 -4a+2o-„,2„,2 , J2 —4a— 4cr„,4\ n /o 

(j + o«(j j [e v i — e ^2 — gy e v i v 2 ' 9y e v 2) = • l^-H) 
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The equation for the inflaton yields 

+ 3«0 + V" - /f (e~ 2a+4 <V + 2e- 2Q - 2< V - 2<? 2 e - 4a+2CT ^ 2 2 - ^e" 40 " 4 ^ 4 ) = 0(2.12) 
From Yang-Mills equations, we have 



and 



f ■ 

v\ + 2y<M + (d + 4ff)u"i + 2c/ 2 e - 2a - 2,T t» 1 t» 2 = (2.13) 

v 2 + 2^^ 2 + (a - 2&)v 2 + g 2 Y e- 2M v\v 2 + # 2 e" 2 "- 2 ^ = . (2.14) 

Note that, when we put v 2 = 0, the above equations reduce to those in abelian cases. 
We need to specify the inflaton potential V(<f>) and the coupling function /(</>) in order 



to solve the above equations. From the previous analysis in |12[ , we know that the condition 

f V 

jy > 2k 2 . (2.15) 

is necessary for anisotropic inflation to commence. When we consider a simple chaotic 
inflation 

V(<f>) = ^m 2 2 , (2.16) 



the simplest choice is [36 



/(0) = e -V/2 _ (2.17) 

Then, the condition Q2.15| ) for anisotropic inflation yields c > 1. It should be stressed 
that the anisotropic inflation occurs for a quite broad class of potential and gauge kinetic 
functions as long as the condition ( p. 15 ) is satisfied p7| 10[ . 

In the next section, we solve Eqs.( |2.9| - |2.14|) and study the anisotropic inflation caused 
by the gauge kinetic function for the SU(2) Yang-Mills field. 
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3 Anisotropic inflation 



Because of the non-linearity of basic equations of motion, it is difficult to obtain ex- 
act solutions. Hence, we first solve equations of motion numerically and find features of 
anisotropic inflation. From the numerical analysis, it tuns out that the nonlinearity of 
gauge fields can be neglected during anisotropic inflation. Thus, we can make an analytic 
treatment of basic equations of motion under the slow roll approximation. 

3.1 Numerical analysis 

Let us solve the Eqs.( |2.l0 
parameters as 

«=1, c = 2, g Y = 0M, m = l(T 5 , (3.1) 

and initial conditions as 

= 12, = 0, V! = 0, vx = 2.47 x 10 -75 , v 2 = 0, a = a = a = . (3.2) 

The initial value for a is determined by the constraint equation ( |2.9| ). We use these param- 
eters for all numerical calculations in this paper since the qualitative result does not change 
even if we change these parameters. However, the initial condition for v 2 , which is not 
included in the above set fl3.2|) , changes the qualitative behaviour of the inflation. Thus, 
in our numerical calculations, we vary the value of v 2 and study features of anisotropic 
inflation for various v 2 , s. 
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Figure 1: Te phase flow in 0-0 space for v 2 /vi = 0.5. We can see two phases of inflations 
which correspond to isotropic and anisotropic inflations. 

In Fig.|], we show a trajectory with the initial condition v 2 /v\ = 0.5 in 0-0 plane. 
Taking look at Fig.|l|, we see that the behaviour is similar to that in an anisotropic inflation 



2.14p numerically. In our numerical calculations, we set 
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in the case of the [/(l)-gauge field |I2"| . There are two phases of inflations, isotropic and 
anisotropic inflations. Since we have started with negligible vector fields, the trajectory 
goes into a conventional isotropic inflation. During this stage, the energy density of the 
vector fields rapidly increases and the inflation soon becomes anisotropic in the second slow 
roll stage. There, the increase of the energy density of the vector fields saturates due to the 
backreaction of the vector field. Therefore, the anisotropic inflation is an attractor solution. 
After a sufficient e-folding, the inflation ends with reheating. 



x 

55 




Figure 2: Time evolution of anisotropy E/if is plotted against e-holding number a. These 
curves correspond to V2/V1 = 0.5, 0.75, 1.33 and 2.0 from top to bottom. We see that the 
anisotropy can be either positive or negative depending on the ratio v^jvy. We also see the 
rapid oscillation of the anisotropy during the reheating. 



In Fig.||, we show time evolution of anisotropy, H/H = a/a, for v 2 /vi = 0.5, 0.75, 1.33 
and 2.0. It is remarkable that the anisotropy can be either positive or negative depending 
on the initial ratio v^j^x- On the other hand, for the anisotropic inflation by a 17(1)- 
gauge field, the anisotropy did not depend on the initial condition for the gauge field and 
it was always positive. In |I8| , it was shown that the statistical anisotropy generated by 
anisotropic inflation with a [/(l)-gauge field has an opposite sign to the one claimed by 
the analysis of WMAP data in the CMB. The negative anisotropy, if true, may suggests 
that the anisotropic inflation assisted by non-abelian gauge field generates the statistical 
anisotropy. We should note that the massive vector fields could also produce the negative 
anisotropy [ £0| . We note that the initial condition dependence does not mean the loss of 



predictability because the initial condition dependence comes into only in the degree of 
the anisotropy and hence it can be absorbed by rescaling the model parameter. Indeed, 
the consistency relations among observables found in are independent on the initial 
conditions. 

In the reheating stage, we see the rapid oscillation of the anisotropy in Fig.|2], which 
implies a rapid oscillation of the gauge field. In Fig.pl, we show the dynamics of the gauge 
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Figure 3: The dynamics of V\(t) and V2(t) with the initial condition V2/V1 = 0.5. The dashed 
and solid curves correspond to Vi(t) and v^t). The horizontal axis is the cosmological time 
t in the unit of k = 1. We can see the chaotic behaviour of the gauge field after inflation. 



field, Vi(t) and ihif) for i'2/vi = 0.5. We can see the chaotic behaviour of the gauge field 
during the reheating. We discuss the Lyapnov exponent which characterize the chaos in 
section||. 



3.2 Slow roll approximation 



In the previous subsection, we found the slow roll inflation with the anisotropy of the 
expansion of the universe. Now, we solve basic equations ( |2.9| - |2.14] ) using the slow roll 
approximation. 

In the inflationary phase, the scalar field takes the value K<p ~ 10. Then, the gauge 
kinetic function has a very large value, f(4>) ~ e 100 . As one can see from the action ( |2.2| ), 
the gyj f(4>) can be regarded as an effective gauge coupling. During the inflation, the 
effective gauge coupling becomes very small gyj '/(</>) ~ e -100 . Therefore, we can neglect 
the gauge coupling during inflation. Then, we can integrate Eqs.( [2.13"D and ( [2.14Q as 



t>i = r 2 (cj>)e 



-a— 4cr 



Pi 



v 2 = r 2 (0)e 



P2 



(3.3) 



where pi and p 2 are constants of integration. Since the energy of the gauge field should be 
subdominant during inflation, we can ignore a in Eqs.( 2^ - 2.10 ). Thus, using the slow roll 
conditions, 2 << V and (f> « V, these equation can be written as 



a 2 



K 



3ad> + m 



-4a 1 



cn 2 (be 



(pi + m)\ 

CK 2 <j> 2 -4a/ „2 



'(Pi 



ipl) 



(3.4) 
(3.5) 
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When the effect of the vector field is comparable with that of the inflaton in ( |3.5| ), namely, 
whenm 2 ~ cK 2 e~ CK2<i>2 - ia {pl+2p 2 ), we find e- CK2<t>2 - ia {p\+2p 2 ) / '(m 2 2 ) ~ l/(c/t 2 2 ) ~ 1(T 2 . 
Thus, we can neglect the second term in the right hand side of Eq.( |3.4j ). Then, from 
Eqs.(|3j) and @, we find 

^ = A + - 2 e- c ^- Aa (pi + 2pl), (3.6) 
da k z m z 

Integrating the above equation, we obtain g-c^ 2 - 4 " — m 2 (c — l)/(c 2 K 2 (p 2 + 2£>2)(1 + 
De" 4 ^ 1 ' )), where D is a constant of integration. This solution rapidly converges to 

CK 2 02 _ 4a= m 2 (c- 1) 

c 2 K 2 (p 2 + 2p 2 ) ' { ] 

From Eq.( p.ll| ), we obtain 3a& = (k 2 /3)e~ CK2 ^ 2 ~' ia {p\ — p 2 )- Therefore, the anisotropy can 
be evaluated as 

£ _ jf_ CK ^_ 4 , _ 2 _ 2(c-l)(p 2 -p 2 ) 1 ,„ 8 , 

H~9a 26 [Pl P2) ~ 3c 2 (p 2 + 2p 2 ) K 2 <f) 2 " 

To obtain the last expression, we have used Eqs.( |3.4j ) and ( |3. 7| ) . From ( |2.9| ) and ( |2.10| ), we 

have ci = — (k 2 /2)</> 2 — (K 2 /3)e _CK2 ^ 2_4Q (p 2 + 2p 2 )- Thus, the slow roll parameter is given 
by 

a 2 . . 

a" 1 cK z (p z 

Therefore, the anisotropy can be written as 

e _ (c-DM-rf) 

i/" 3c( P ; + 2pI) 1 ' 



We should notice that the abelian result can be recovered if we put p 2 = |12| . From the 
expression ( |3.10| ), we obtain an inequality as 

C_1 6<^<^6. (3.11) 



6c ~ H ~ 3c 

It implies that H/H is suppressed by the slow roll parameter. The anisotropy can be either 
positive or negative depending on the ratio p 2 /pi- In particular, when p\ = p 2 , we have no 
anisotropy. Except for this accidental case, we have the anisotropy proportional to the slow 
roll parameter. In particular, for the cases p 2 > Pi, we have a negative anisotropy H/H 
which never occurs in abelian models. This analytical result explains the numerical result 
shown in Fig.|2|. 

The statistical anisotropy induced by an anisotropic inflation can be characterized by 
the direction dependent power spectrum 

P(k)=P (A;)[l + ^(k-n) 2 ] , (3.12) 
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where Po(k) is the isotropic part of the power spectrum, k is a wavenumber vector of 
fluctuations and n is a specific direction, in our case, this is the x-direction. Here, the 
number g* is a magnitude of the statistical anisotropy. According to the recent analysis of 
the CMB data JTTJ . (/* takes a positive value. The prediction of anisotropic inflation based 
on the abelian gauge field models was inconsistent with this result []15-|18||. In our non- 
abelian models, however, the anisotropy can give rise to a positive g* since the anisotropy 
of the expansion can take any signature depending on the ratio "Pil"P\- Thus, if the CMB 
data shows the statistical anisotropy with positive g*, it may imply the anisotropic inflation 
with a gauge kinetic function for a non-abelian gauge field. Of course, we should keep it in 
mind that the data analysis may contain systematic errors. 



4 Chaos during reheating 

In this section, we study the chaos during reheating in FigfJ. The chaos in Yang-Mills 



cosmology has been studied previously [42-44]. In those analysis, the existence of coherent 



non-abelian gauge fields is assumed. Here, the initial coherent non-abelian gauge field is 
provided by the anisotropic inflation. 



As we explained at the beginning of subsection. |3.2j , the gauge coupling is extremely weak 
during the inflation. Hence, we can neglect the non-linearity of the gauge field. During this 
phase, the coherent gauge field is produced due to the rapid variation of the gauge kinetic 
function. While, at the end of the inflation, the gauge kinetic function becomes f(4>) ~ 1 
and effective gauge coupling gY/f{4>) becomes of the order of gy. The chaotic behaviour 
after the inflation occurs due to this relatively large effective gauge coupling. Therefore, 
without the time varying gauge kinetic function, this chaotic behavior of the gauge field 
never happens. 

From the numerical analysis, the pattern of oscillation in Fig.[3] seems to change depend- 
ing on the ratio Vijv\. The anisotropy H/(eH) also depends on the same ratio. Hence, we 
might have some relation between the anisotropy and the chaos. So, we need to check if the 
relation between anisotropy and chaotic behaviour exists. We consider linear perturbations 
of Eqs.( ^9| - p.l4D by the substitution a — > a + 5a, a — > a + 5a, (f> — > + 5(f), V\ — > v± + 5v\ 



and V2 — > i>2 + 5v2- Since the chaotic system is sensitive to initial conditions, the per- 
turbation would grow at late time. We define "Lyapunov exponent" X(t) as 5v2 oc e A ^*.[] 
The X(t) represents the growing rate of the perturbation. In Fig.[|, we depict the X(t) for 
£/(e#) = 6.20 x 10~ 2 , 2.26 x 10~ 3 and -7.10 x 10~ 2 , respectively. We cannot find any 
correlation between X(t) and S/(eiJ). It indicates that Lyapunov exponent does not de- 
pend on the anisotropy in the anisotropic inflation. Although there is a gauge ambiguity 



1 Mathematically, the Lyapunov exponent A is defined by A = limt-»oo ln(|iw(f)|/|io(0)|)/i where 
w(t) = (6a,- ■■ , 6v2,Sa, ■■■ ,5v2)- Since it is difficult to solve Eqs.(2li|- 2.14 ) for sufficiently long time 
in our numerical calculations, we define the X(t) in this way. 
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Figure 4: The Lyapunov exponents X(t) for various initial conditions. Solid, dashed 
and dotted curves correspond to E/(eH) = 6.20 x 10" 2 , 2.26 x 10~ 3 and -7.10 x 10" 2 , 
respectively. We cannot find any correlation between the Lyapunov exponents. 



in defining the Lyapunov exponent, our conclusion itself does not depend on the choice of 
the gauge. 



5 Conclusion 

We have studied the anisotropic inflation model inspired by supergravity where the 
gauge kinetic function for the SU(2) Yang-Mills field is non-trivial. We found that the 
anisotropy of expansion rate can take either positive or negative value depending on the 
initial ratio v%jv\. Namely, the shape of the comoving volume becomes either prolate or 
oblate depending on the initial configurations of the gauge field. This new feature can be 
attributed to the multi-component nature of non-abelian gauge fields. In principle, this 
could occur even in multi-abelian models. However, it would be difficult to organize the 
multi-abelian fields so that the same result as the non-abelian models can be obtained. On 
the other hand, the gauge structure of the non-abelian gauge field self-organizes the config- 
uration. In spite of the above initial configuration dependence, the anisotropic inflation is 
still an attractor in the sense that the details of initial conditions are irrelevant except for 
one relevant parameter which controls the anisotropy of the universe. In addition to the 
anisotropic expansion, we found the chaotic behaviour of the gauge field during reheating. 
This is due to the nonlinear self-coupling of non-abelian gauge fields. We calculated the 
Lyapunov exponent of the chaos and found that the Lyapunov exponent does not correlate 
with the anisotropy of the inflation. This indicates the universality of the chaotic behaviour 
of Yang-Mills field in an anisotropic inflationary scenario with a gauge kinetic function for 
a non-abelian gauge field. 
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It is remarkable that the anisotropy E/i? can be negative in our inflation model. In ||18|| , 
it was shown that statistical anisotropy generated by anisotropic inflation by [/(l)-gauge 
field takes an opposite sign to the one of observationally favored value. While, the negative 
anisotropy would lead to the observationally favored signature. It may indicate that the 
anisotropic inflation with a gauge kinetic function for a non-abelian gauge field generates 
the statistical anisotropy in the present CMB. Probably, it would be too early to conclude 



something, at least, we should wait for more precise CMB data provided by PLANCK [45 



It is also intriguing to observe that the new relevant parameter v 2 / V\ helps to weaken the 



constraint on the model parameter c found in ||18|| . This is because we have two parameters 
v%jv\ and c to control the anisotropy. 

In this paper, we have assumed axial symmetry for simplicity. It is a straightforward 
exercise to consider more general configurations. As a future work, we can consider non- 
gaussianity in an anisotropic inflation along the previous paper p0fl . From this point of 
view, it is interesting to extend anisotropic inflationary scenario to a non-slow roll type 
inflationary scenario such as the DBI inflation model fHfl . It is also interesting to study 
implications of chaos in the early universe. For example, gravitational waves may be gen- 
erated at a reheating phase. Because of the chaotic behaviour of the gauge field, we may 
be able to find relic of the chaos during the reheating. 
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A Axi-symmetric gauge fields 
A.l Symmetry constraints 



Using symmetry, we can make a reduction of variables [f0[]. In the present system, 
there are spacetime isometry and gauge symmetry. Using these symmetry, we will make 
the variables as simple as possible. 

Imposing the translation invariance along d x , d y and d z , we can put the gauge field as 
A = A t {t)dt + Ai(t)dx l . Furthermore, using the local SU(2) gauge freedom, we can fix the 
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time component of the gauge field as A t (t) = 0. Then, the gauge field can be written as 

A a = f3 a (t)dx + Y{t)dy + 5 a (t)dz . (A.l) 
The residual global gauge transformation is given by 



i[A,u] a = e abc u b A c 



Mbc„.b 



u\P c (t)dx + i c (t)dy + 5 c (t)dz) 



(A.2) 



= [u x (5) a dx +(ux ^) a dy + (u x 5) a dz . 

where w a are constants and we used a vector notation like as u = (it 1 , it 2 , it 3 ), /3 = 
(/3 1 , /3 2 , /3 3 ), etc. Without loss of generality, choosing an appropriate basis in Lie algebra, 
we can put u = (u 1 , 0, 0) = u l e\. 

Now, we consider the rotational symmetry. The rotational transformation is generated 
by Crpdx = 0, Cc/jdy = —dz and C^dz = dy. Under the infinitesimal rotational transforma- 
tion, the gauge field transforms as 

C^A" = -7 a (t)dz + 5 a {t)dy . (A.3) 

For the rotational invariance, the above rotational transformation must be absorbed by the 
gauge transformation (|A.2j), namely, 



L^A a = 5 g A a . 



(A.4) 



Since the gauge field has to be the same after making a round, A = exp^nC^A holds. 
Thus, using Eq.( |A.4j ) and [C^, S g ] = 0, we derive a relation 

A = exp(27rC^)A = exp(2ir5 g )A . (A.5) 

Therefore, exp(27r5 9 ) must be an identical transformation. This implies u 1 = n £ Z . Then, 
substituting Eqs.( |A.2j ) and ( |A.3| ) into Eq.( |A.4j ), we obtain 

(A.6) 



For n = 0, we find 



ne[ x f3 = , 
ne[ x 7 = 5 , 
ne[ x 5 = —7 



/3 = (/3 1 (t),/3 2 (t),/3 3 (t)) , 7 = 5 = 0. 



(A.7) 
(A.8) 

(A.9) 



For n = l, from Eq. (|A.6 ), we obtain /3 || e\. From Eqs.( [A.7|) and ( A.8|) , we find e*i _L 7 _L 
5 1 e*i. Furthermore, taking the absolute value in Eq.( |A.7 ), we have |t*| = Thus, the 
vectors /3, 7 and 5 can be written as 

/3 = (/3 1 (t),0,0) , 7 = (0, 7 2 (t),7 3 W) , 5=(0,- 7 3 (t),7 2 W) • (A.10) 

For n > 1, Eqs. (|A.7|) and ( |A.8|) cannot be satisfied unless 7 = 5 = 0. Thus, we obtain 

/3 = (/3 1 (t),0,0) , 7 = 5 = 0. (A.ll) 
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A. 2 Yang-Mills constraints 



In the previous subsection, we classified the gauge field into the three types (|A.9 ), (]A.1Q ) 
and ( A.ll ). Here, we impose the Yang- Mills constraints on these expressions. 

First, we consider the case ( |A.9|) . Substituting Eq. (|A.9|) into a time component of 



Eq.([2.6|), we obtain 



?2nl 



?2 o3 



3/o2 



j3ol 



?1 fl3 







(A.12) 



Therefore, we have /3 = /?i(t)c, where c is a constant vector. Thus, choosing an appropriate 
basis in Lie algebra, we can put 

/3 = (/3 1 (t),0,0) , 7 = ^ = 0. (A.13) 

Next, we consider the case (|A.10|) . Then, from a time component of Eq.( [2.6|) , we find 

7 2 7 3 -7 3 7 2 = 0. (A. 14) 

Thus, we get 7 2 = C7 3 , where c is a constant. Choosing an appropriate basis in the subspace 
of Lie algebra spanned by T 2 and T 3 , the vectors /3, 7 and 5 can be written as 

^=(/3 1 (t),0,0) , 7 = (0, 7 2 (t),0) , 5=(0,0, 7 2 (t)) • (A.15) 

For the case ( |A.11| ), constraint equations in Eq.( p.6| ) are trivially satisfied and we have the 
same expression as Eq. flA.l3p . Therefore, Eq. ( |A.l5| ) is the most general expression for an 
axially symmetric gauge field. 
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